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REVIATIONS AND SIPEOLS TO FE USEN IN TiS: REPORT 


Wing aspect ratio 

Aircraft total draz toef/icient 

Induced drag coefficicnt E 
Drag coefficient at zero lift 

Thrust specific fuel consumption 


. Centrifugal force 
Total aircraft drag 
Total differential operator 
Partial differential operator 


Induced drag 
Drag at zero lift 


Airplane efficiency factor 
Total effective thrust 


Thrust due to momentun change of air 


Effective thrust available 


Gravitational acceleration constant (32.17 ft/sec?) 


Altitude 

Aircraft lift ` 

Mass of air and fuel 

Mach numler 

Dynamic pressure (q = 1 py%) 

wing surface area 

Time 

Engine exhaust velocity relative to the origin 


Exhaust velocity relative to the Lody (V - u) 








= 
Velocity 

Speed of sound 

Velocity; derivative with respect to time 
Total weight of Lod, or aircraft 
weight of fuel expanded from point 1 to 2 
dorizontal distance 

Flight path angle relative to the x-axis 
Flignt path an:le derivative with respect to time 
Air density 

Angular velocity in vitch 


Aircraft total lift coefficient 
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and altitude, such as nt. A on fig. 1, and flown a [piven distance X]. 
Tue schedule of altitude and sneed would Le detorrnined ty the inter- 
relations of the wrialles of the proLlem bounded Ly the capalilities of 
tne particuler aircraft. The problem coule just as well ke started 
forn the instant of takeeoff, say Dt. Ps and thus co the entire flight 
fron the takesoff until over xj. 

Commencing a maximun r: nge or minimum fuel schetul> from 
take-off is a completely new concept. Classical renge solutions do 


not attempt to solve for an optirium flight path from point of take-off 


a e e * [] af” 
to optimum cruising altitude, Lut ratner specif;, for turbo-jets, 


P 


amnebulous "minimum fuel clim" to tne highest altitude ət which 
speed for (Cx'*/CD) max is equal to Vj. The aircraft then cruises 
at taat altitude, slowly climbing as the weicht is decreased, thus 
increasin: tno altitude at which the speed for maximum ci Jon is 
equal to the maximum speed. 

This "minimun fuel climb" will, if properli deterninec, 
achieve a certain aitituds with a minirum fuel exenditurs, Lut this 
is not the primary concern. The proller is to travel a distance x} 
with minimum fuel consumed. If the variables of the problem call for 
a climL to a specified altitude, then that climb should be executed to 
give maximum distance for a viven fuel expenditure to get to that 
altitude. As long as the fli, ht »»sth stars within the ¿imitations of 
tno aircraft, however, there snouid be a single path with no dis- 
continuities for ninimun fuel path for a given distances, Thus 


looking at fig. 2, it can be seen tat to travel from pt. a to some 


pte ¢ 2 distance x] away there must be some minimum path ac which uses 


N 











less fuel than the path abc which is composed of two parts, even 
though ab ani Le are separate paths. 

Classical range equations, moreover, are not generally 
applicalle or accurate for present Jay high performance aircraft. 
Classical range performance is based on assumins that the flight 
path angle (1) is small such that its cosine is considered one, 

(this assumes that L = W) and that flight path accelerations along 
the path were immored (i.e. Thrust = rag). From a review of the 
literature on the sulject, it was determined that everything written 
to date on the rante of aircraft nes been based on these simplifsing 
but very restrictive assumptions. With modern hizh performance eir- 
craft, flicht path angles un to and including 90 degrees are vossille 
and flicht path accelerations can easily exceed 1.0 Ze In fact, in 
considering the »ener»l probiem of pouere? flicht (rockets, missels, 
WO, etc.), the convenient and conventional condition of steady lel 
flisht is almost the exceotion instead of the rule, 

It is the intent of this work to indicate methods based on 
the principals of variational calculus for tne det:rmninetion of schedules 
for minimum fuel or maximum range flights for all types of rowered 
Tlithts. It is intended, more specifically, to develop equations 
Bee vie minimun fuel paths of turbo-jet type aircraft. This differs 
ROM the. classic problem in vs the scheduie Lerins at teke-ofi and 
minirizes the fuel expenditure ror that point to a distance xi, taking 
ints consideration possible large values of flight path angle and 
consicershle flisht path accelerrtions. In other words, the problem 


consiccred nere is the broad ıreneral proLlem. It uses accurate, 


realistic conditions ani ninirizes fuel consumption for the entire 








flit path from "wheels up" to "over target". 

It is hardls necessary: to discuss the inmortance of reduced 
fuel consumption since its advanteres are so readily apparent, but 
to help keep these advantages in mind, a few of then will be listed: 

(1) Increased useful paload (bombs, cargo, passengers, etc) 

(2) Creater range or combat radius with the same payload 

(3) More economical operation overall 

(L) Complete specification of a flight schedule which could 

be incorporated into an automatic control for a manned 

or pilotless aircraft. 
The adeptability of these and other advantages to all phases of 
aviation is also readily apparent. The long range bomber, ‘the short 
ranze fighter, the commercial carco and passenger carrier, the pilotless 
aircraft, and practically all other vehicles traveling through the 
air could adopt flight schedules determined in this manner to increase 
their effectiveness and their efficiency. 

A specific application for which this work could be im- 
mediately useful is tne determination of entry and escape paths for a 
supersonic bombing missione. Those parts of the entry and escape 
patà not determined b; other tactical considerations wouid most prob- 
ably be flown at fuil power including afterburner and following a 
path to give best fuel economy at maximum power. The bounded thrust 
path equations ceveloped laterin this thesis apply directly to this 
probleme 

In this and many other applications, Mach effects must be 
considered. For simplicity, however, Mach effects have been omitted 
from the derivations that follow, but have been treated separately in 


Appendix Be 
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PART TWO: DEVLOR ENT OF “SQ ATIONS 


The range of an aircraft is, as we have hypothesised, a 
function of its flicht path. Its flight path, in turn, depends on 
variables such as amount of fuel, velocity alone path, fliyit path 
ande, accelerations along and normal to the path, and thrust spec- 
ific fuel consumption. These varialles are, in turn, functions of 
each other, drac, thrust, lift, weight, and air density. Relating 
these variables among sufficient equations to solve for a minimum 


fuel path is the problem of this thesis. 


A) EQUATIONS OF MOTION 

The equations of motion pertinent to the problem are the 
force balance equations of a body moving in a vertical plane. This 
motion will include accelerations- along and normal to the flight path, 
as in maneuvering flisht, since these accelerations will, a priori, 
affect the loads and thrust required, and therefore affect the range 
(or fuel consumption). The forces acting on a body moving along an 
arbitrary path in two dimensional space are shown in fig. 3. 

If this body moving ip space is a conventional aircraft, 
the forces Fy, Fo, and F3 can be expressed ası 

F4 The total effective thrust. (F) 

Fo The lift of the moving body. (L) 

F3 The total drag force. (Body drag, D, and inertia, m E ) 
Ey equating forces parallel and perpendicular to the flight sath to 
zero, the followiny equations are oltained: 

(1) Forces parallel to flight path 


DE “OeF1 -F3-Wsind (1) 
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CF 





or for an aircrafts: 


F*-D - i: mV) - Wsin3 + u — = 0 (2) 


where F, D, W, 3 p ü, and m were defined previously and: 


d(mV) = m dV * V dn 


dt dt dt 

But 
mur 3) 
dt ws di 


EF, = F®-p - mV + VCE -~ Wsin3 - uz = 0 (h) 
g 8 


Considering the definition of thrust, it can be Seen that the total 


effective thrust acting on the aircraft would be equal to that im 
parted by the air through the engine as it experiences a change of 
momentum plus that imparted by the fuel mass leaving in the exhaust 
due to its momentum change. The thrust due to the leaving fuel is a 
function of the difference between aircraft velocity and exhaust 


velocity and is.expressed by: 


(u =- V) Cp F | 
E 
Thus ! 
F = F* . (u =- V) CF 
8 
= F* + ue Cp F (5) 
8 





The equation now Lecomes 


F- D.N V- hsini -0 (6) 
E 


(2) Forces perpendicular to the flight path 


LF, s (0 = Po - CF - W cost (7) 


Í 
or for en aircraft 


L- ve =G (8) 
But 
CF =m VW 
and y 
us - di = Y 
dt 
then i 
CF =wWVY (9) 
g 
And therefore 
PRo-L-WVÉ-Weof -0 (10) 
g 


e 


Equations (6) and (10) above are the general equations of 


motion pertinent to flight in the vertical plane, These will be modi- 


fied later for use in the range problem to the extent that derivatives 
with respect to time will be replaced by parametric derivatives with 


respect to either weight or distance. 


B) EQUATIONS OF MOTION FOR RANGE 
Before attempting to maximize the general equations of motion 
we will further modify them to be more applicable to the problem of, 


(1) Maximizing the range, x, for a given fuel expenditure, 1W2s or 


(2) Minimizing the fuel expenditure, jWofor a given range, % 


n 










Be fe of Saye 
v di “at dt 


e 


wu uw 
Substitutino (11) and (12) into (6) an' (10) rives: 
L suv gt Cp * - eos} = 0 


GQ 


P-DéMVY op F - Ysind = 0 
R 


The integral equation relating x andi in simplest fora is: 
xs dx d 

d 
ge = dx gt 
A dT 
a 


But = Veos$ 


dr CT i 
therefore 
Zw ~ Vcos di 
Y 


any 
pe (13b) ER 
= 


Po. ¡Due "sin = cos! - 1 
f 1 + CHO, Er y Cr 
g 8 


Tor maximizing ranse, X, it is convenient to chanse variables 


(11) 


(12) 


(13a) 


(13b) 


(16) 
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10 








Substituting into (15) 


x 2 y Vo y! cos Y dw 


E (17) 
L --Weos ¥ 
or 
x = - VeosY (1+ Vy'op) du 8 
1 
Cp (Dwsin (18) 
(2) For minimizing the fucl used, jko, the changes of 
variables found useful are: 
Definine a( ) $ ( )' 
Then 3 
V = dV = dV dx t ViVcosY (19) 
d "dx at 
And l 
$= d = dd dx = X'vcosd (20) 
© metat 
Substituting into (6) and (10) gives: 
L-W ve Y! cos =- Wcosf- O (21a) 
B 
t € 7 = 
Es oe V! Vcosý =- Wsin¥ = 0 (21b) 


The integral equation relating x and W in simplest form for minimizing 


1"2 Js: 
W= dw 
m dx 
As before 
dil > CT F Vet (22) 
dx Vcost 
and 





= 
N 
"lm 
< 
olv 
O 
ui zj 
A, 
E 








1? 


Solving (21b) for F and substituting gives 


a or (D + + Wy'Vcos¥ + WsinY)dx 
V cos 


we |. D+ WV + Mal 
J u = me | ax (23) 


Now, although equations (18) and (23) relate the variables 


or 


of the problem, the: do not specify the flisht path. They are equations 
for the general motion of a body in the vertical nlane, expressed in 
such a way as to give the relationship between range and fuel ex- 
penditure for any flight path. The variables under the integrals 

are dependent on the path, yet the equations as such do not specify 

a path. To determine equations for the particular path which will 
give the greatest range or least fuel expenditure it is necessary 

to look to the mathematics of maxima ani xd Presenting the 
problem in more mathematical terms, the unique path desired is that 
patn giving the extremal value cf an integral. The problem of finding 
an estremal value belongs to that branch of mathematics known as "The 


Calculus of Variations". A brief discussion of this subject follows, 


C) THE EULER EQUATIONS OF THE CALCULUS OF VARIATIONS 
Let I be a quantity defined as an integral function of 


Several variables as follows: 


mo sen, ---,3',2!, ---,) dx 
et 
where ( )' indicates the first derivative with respect to x, The 
Calculus of Variations is concerned with finding the extremal values 


of the quantity so defined, The necessary and sufficient conditions 





ee te — A EE — —  —— n= a a 


Fig A 


FAMILY OF PATHS FROM 
POINT ON BON 





Yo (2) = Y (x) + € py" (x) (26) 
With the conditions: 


Tas) * TH 


N i.e. 1 (x) = O at the end points 
Y 0) 7 ¥(x}) 


Where €, is a small constant quantity and no?) is an arlitrary 
function of x that is continuous and has continuous first and second 


derivatives in the range of integration, 


Y, (x) therefore defines all possible paths of integration 


between X5,Y4 and Xj,Y, so that (2h) may be written: 


XY, E. u 
In * [° [x5 Yn G0, Ys! (x) dx 
Y, 


0) 0 


or replacing u and Yn’ from (25) and (26) 
X,Y 
I,* [alx,v(x) EY D PE N'(x)) dx (27) 
Xo, X, 
It is.clearly obvious that Ig differs from Ik only through the value 
of € and it may be concluded that: 
Ip“ I€) = f(e) only 


and since € is not a function of x, the ordinary maxima ahd minima 


rule may be applied, 1.e€¢3 


d =Q for an extremum 
dE lero 
and, furthermore, the differentiation may Le made under the integral, thus: 


Y, Ñ 
d = da ‘ | 
de i de dx - 
X 


®©) % 





From (27) 


dG = AG e d| ¥(x) + enví(x) 


de dr) + En (x) de 





+ AG e dly' (x) +16 y no $0 di 
oly (x) + €, y" (x) | de + dé 


Since 


dx = dy 2 dy = dY'! = dy! = Q 
y >" Y 
then 
m open 7 (x) 
3 [Go + G2] JE SË + on (x)] 
KY) 
Or 
di - y 
de t (et T Ae (29) 
Xor Yo 
and 
XY X, Y, 
mt O mu Bean. ng 
de o(Y * En?) AP + as ") 
Ko, Yo Xo, Yo 


From whence 


m 
oG (ex -- 30 Vi 
E OT, Një 
Integrating = right side E =. and remembering "that 7 E.A) 


7 en» n) F dr m 


' - a[2G 
OG. Idx Së (a) (30) 


dyn! " dx 
X, fe X, Yo 
Substituting (30) into (29) gives: 
X,Y, 
di E JË + dd 396 3 (31) 
Xor Vo e 


Evaluating (31) for E= 0 to find the extremum: 


n 
òG -d 36 $ 
fts & dia SË 


X, Y. 








But this is true only if 


i8 aer 


throughout the range of integration, Since Y is completely arlitrary 
at every point except the end points, it necessarily follows that for 


an extremum: 

¿($ e (32) 
Y daxldY 
This equation is called the Euler equation. It must be satisfied if 
the integral is to have an extremum value and it is through this 
satisfaction that the path, Y(X), is found. | 


Proceeding to the more general problem of several varial les, 


the integral must satisgj the definitions as before and: 
: XII, 


I= /G|X,¥,2, ---, ¥,2', ----] de 


Bo ^" o 


Proceeding as before, we arrive by an analagous developement at the 


Euler equations: 


$3 m | (33) 


These equations (33) must be satisfied sinultaneously for an extremal 
value of the integral, such satisfaction therefore defining the desired 
path as in the two dimensional case, 

It should be remembered here that an extremal value for the 


integral may be either a maximum or a minimum value. Thus in deter- 
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dr 


mining a solution for equations (33), we must determine ly other E 4 
means whether a particular solution will yield a maximum path or a 
minimum path. For most types of enzineering work this will be in- 


tuitively apparente 


D) THE EULER EQUATIONS FOR MAXIMUM RANGE 
86) Minimum Fuel Zquations (i.2., x independent variable) 


Equation (22) is of the fom 


X2 
- WW = I(c) ‘|= 
të 


where I(c) = -W the quantity (of fuel) to be minimized and G = Cp F 
Veost 
the integrand function, 


Now: 
G = G(C7,V,2,F) 
F » F(D,V! ,u,V,d) 
D = D(L,V) 
L = L(W,V,3,81) 
thus 
G = G(W,V,3,0p,8',V!) 


j 


Now, Crs the tnrust specific fuel consumption, is a variable, and a very 


important variable where range is concerned, but treating it as such 
in the derivation of Zuler equations gives rise to extremely com- 


plicated relations and to derivatives on Cp that are very difficult 
vo evaluate. We have therefore elected to treat G as a stepwise 


variable constant, 


Thus G = G(W,V,3,V',8') 


19 


And the Euler equations of this function become: 


ZEE 


20) 29" 2 (87) = 


3m 
3.) 00 - d/oc z 
m al) 


le) Derivation of the first Euler 


- (ge 
dx \ dV", 
Now from (22) 
G = O4 F 
Veost 
Then 
4 s _Veosé Op Sy - Op F cos% 
(Vcos$)? 
=G E E aj 
F V 
And 


ðG = CT OF 
OV! — Vcos? ES 

d òG = Vcosi Or OF 9F ` V'cosY - Vd sing) 
ax ove le g ove IM - 


* (Vcosd)2 


Using an abbreviated notation: oi Y 
dié 2 


AG 2$) G Fyr = F y: (V! cos Y - Vg' sinti) 
dx VOV Fo F Vcos % 


or 
Gys” = G| Fy! -= VER, + = Pim] 
3 vr. 
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-GFy- |Py:' = Vi + Pw 
Tin v : 


Now 1 


/ 


Gy = Gyv! ~O =G)Fy 1] -GFyt EN - Y! S tan) 
F v F Fui V | 


Multiplying through by F gives us the first £uler 
t G 


Fy -F =- Fyı! 4 VE Fyı - Fyr X'tanë t 0 (3h) 
Y y 


2.) Derivation of the Second Euler equation 


dG- d dG = O 
ax 


as befores 


G = Cp F 
Ycos $ 


OG = Vcosf ers + Or F Vsind 
“38 (Vcos®d >) 


and 


then 
d/ac a Fr ' - Visi 
aS) = C Vcos? I Fy (V! cos $ VÀ sinf) 


 (veoshë 0 


BE E ja - Fy x u^" 


= G Ein - Vv! sv 










que F 
- 
g the Buier: 


f Mi za de * *tany]- G Py PRA -de qty 
IF + T v 


Multiplying through by F/G leaves the second equation: 
Fy + Ftant - Ps! 9) V' Fy: = Fy Y 'tany - 0 (35) 
3.) Derivation of Third Euler 


ds uo 
OW drit 


G = Cr F 
Ucos Y 


oW NË) 


dG = ¿Ju!=1 j 
dw dur 


d (që » O 

dx (dit 

Gy -~ Gwi’ = G| Pw = 0 (36) 
> 

Therefore Fy > O and there is no Third Euler in this case. 


b) Maximum Range Equations (i.e. Wis independent variable) 


Notes: + ) = ( y! i 


From equation (15) 


Ge. VeosY 
GJË 


1.) Derivation of the First Euler Equation 


H Ae ON yw 
a ae : 
= ern de 
V F | 
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d /dG \= pe (v YE, * Veos TER! VY! 
a (Sir) em qi * opas te 
~/Vcos¥ Fy: 2F F' Cp m 
(Cp F2)2 


= VI Fui cost * Vecos% Fy! - Fui V 3' sink 
Cp Fe Gr F2 Cp Fe 


- 2Vcos¥ Fy Ft 
mr: e 
Cp F 








Quit = -G P Fy + Fyt = ¿itang - =| 
VE me p e om 


: geht e -4 


> F Fys 


Forming the Euler 


Gye Gy = all - Fr |- GFy ja + Vitanf - Fu? - 1 |=0 
p E + AE F Xp 


And 


t 
Fo Fy e CF vy Fi = F yı f! an * Fy! + Y En = 0 I (31) 
Em au. NV 








2.) Derivation of the Second Euler Zquation 


dG - a /dG\= 
sy ° a Ge)" ° 
9G = _V [-F cin’ - cosyY (dF/dY ) | 
PIS Op ye 
= wl} tany + Fy 
d 


AG » VcosY |. SE/AY ' = Vy cos Y 
"ar Cr | (F)2 | ae 


d 28.) 7 - 1 [F,cosY V * Vcos d F,// = VF, Y'siny 
— [YE 41% X à 
cw Vox! e) F2 Um F2 


- VE, cos Y (2F Ft) 
Fl 


n. TË = aha -<G [1 Ey + Fy! — që Fytanë m^ 








F VF F F F2 
Now 
Gy e Gyt = -S|tany* Fy | * G| V Fg, + Fy) - YE tang 
FD (rr rë 
- a = 0 
Fe Ja 
‘From which: 
F tanë t Fy- y Fy - Fp! t PF tani e 2F FI = 0 (38) 
TE 


Thus equations (3h), (35), (36), (37), and (38) represent extremal 
paths in terms of flight path angle (4), velocity (V), horizontal 
distance (x) and pounds of fuel (W). The first three equations re- 


present extremal paths for minimum (or maximum) fuel consumption for 













reling a given distance x, and the last two equations represent an 
extremal path for maximum (or minimum) distance x for a given fuel 


expenditure We 








PART III 


THE POUNDED CASS 


The equations derived in section D of Part II apoly to a 
more or less general propelled body with unlimited thrust and thus 
unlimited velocity and flight path angle. They do not consider any 
boundaries on the thrust available to follow the extremal path that 
the. desimnate. This extremal path maj very well call for altitudes, 
velocities, and accelerations above the physical ability of a parti- 
cular airframe - engine combination for all or part Of the path. © 
Taosg Reef tions are, therefore, ee extremums for an ideal 
body without limitations imposed on any of its parameters, They 
define a path that is the best possible and should be followed when- 
ever that path is within the operating capabilites of the particular 
propelied Lody. 

The primary limitation imposed by physical factors in any 
particular case is that of thrust. Thrust varies from zero to a 
finite maximum value that is a function of the individual thrust 
EE. and other variables. For a given.type power plant the 
maximum effective thrust is usually a function of the variables of 
the problem, For example: The thrust of a rocket is a constant; the 
thrust of a ram jet varies with velocity and altitude; the thrust of a 
turbojet varies with velocity and altitude, Therefore, for any riven 
airframe - engine conlination there is a limiting envelope or Loundary 
of thrust which limits the path variables such as velocity, flicht 


path angle, altitude, etc. Whenever the general path equations of 
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Bert m sexi a thrust that exceeds this boundary, the path becomes 
unrealistic and of antic interest only. It then becomes necessary 
to compute a new path from this — that will be an extremal based 
upon a maximum continuous value of thrust available. Euler equations 
will now be derived to define an extremal beth ter the case of Lounded 


tnrust, 


A.) The Bounded Case for Minimum Fuel for Turlojets 

As was Said before, the thrust of a turbojet varies with 
altitude and velocity, i.e.: 

Be f(V,h) where f pertains to bounded thrust 


But since altitude is a function of path angle, we may also 


say? 
F = £(V,$) 
From before: 
G * OT F 
eos% 
But now 
G = OT £ 
Vcos $ (39) 


This function is independent of the characteristics of the airframe - 


engine combination and must be modified to include them, From (21b): 


V cosý = g(f - D- Wsind) 


WV?! 
Substituting this into (39) above 
C = OfWV (Lo) 


2(f=D-Wsind) 








lo) Derivation of the First Euler Equation 


T 


“6 al IL/AV_ - àD/dV | (là) 
f - D-Wsin 
6 Op f W 


| vt g(f - D - WsinY ) vt 


| But 


ee q 


| G = Wi 
Then 
6 = wt 
“SYT VI 
0 
" T (12) 
da paa = VWa WW 
EE) (y) 
Forming the Euler (11) - (12) = O 
\ | 
a oG - d (9G | - W'[fv- fy - Dy ~ yi _ wem 0 
j Ov dx ( = E m (v1 )? 


| From which 


We Vify =- fy- Dy -jan wv s0 : (43) 
UN f f-D-Wsin WV] 


| 2.) Derivation of Second Euler Equation 


k "la tX - - ànA* - — (uh) 


a - D - Wsin 


ae 
ot 
2G - 


(£-D-WsinY) 


“|e 


SY (f-D-Wsiny)e (£-D-Wsin3 


| t "| Dv | (45) 
| - «Sin 


0% 
2 26) = Ww! EE - Dy (f'-D'-wW'sing - Wy! | 
d 
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Forming the Euler (Lh) - (15) = 0 


y Dy' ë Ly (E - Dy - Wcosg) 
(T-5-V sint) f (Fo DASHI 


- 


4 
| "v. Dy” - Dy (FI Dry! sin - Në cos Y) ate 
Së 
: (f-»-xXsin?)4 
From vhicn 
w Dy! = fy (f-DeWsind) + fy - Dy = Weos¥ + Dy 


TI f 


- Dy LD sing -WÀ! cos Y D 


(I-D-WslnA) 16) 


Partial derivatives of drag with respect to V' and thrust with 
respect to $' are zero since drag and thrust are not functions of 


V' ani $' resoectively. 


E) The Pounded Case for Maximum Range (i.e. with W as independent 


2 





variable) 
As before 
G = Veos% 
| z 
and 


F = f(V,h) = f(V,j) only 
so 


G = Vecos% ) 
"T- eu 


Introducinz the airframe charactcristeics: 








From (13a) 


L¢#Wv$'Cr F - Wcosb =0 
8 
cost = L + vil of 
B 
Substituting this into (7) leavess 


a 


G = VL + pag EX (18) 
x Cpu f g 


1.) Derivation of First Euler Equation 


B = L + Vle . Vint + 2 VY' (49) 
QV Cf CE CW f2 g 


OG = _Viye (thrust is not a function of V') 


E 


d i») 2 ULyr + VI Lys ONES Ly: (Wf! + f) (<0) 
dW \ gir “QUË Gy (uf)2 


Forming the Euler (19) - (50) = 0 


L + Vy = VL fy * 20; 8B £V - Vigo! - Vilyy + Wy [wer + fjeo 
= E wt 
(51) 


2.) Second Euler Equztion 


MG = f m . 18 j 
alte] i E 
òG & VL y! + ve 

au Cyilf E b 








Porming the Euler (52) = (53) = 0 


f 


f 


" | 
Ly -L£y > Ll - Y Ly + Ly[ne! +4] - 200.2 v! = 0 


8 





cH — ES 010 
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PART IV: EVALUATION OF THE DERIVATI VES 
AND FORMULATION OF THE FINAL EQUATIONS 


The Euler equations as developed in Parts II and III contain 
derivatives on several of the variables. To determine the desired 
path, these derivatives must be evaluated and included or discarded 
depending upon their effect. In this thesis only the derivatives | 
for the minimum fuel path, both Lounded and unbounded, will be deter- 
mined. The evaluation of the derivatives for the maximum range path 
(Ww independant variable) is analogous but the expressions are much more 
unwieldy. p 

Considering first the derivatives to be evaluated for the first 


Euler in the unbounded case, we haves (1)Fy, (Fyr (3)Fyı!. From 


. equation (21b) 


F=D+wv — WsinY (55) 
g 


Taking the desired derivatives: 


Fy =Dy + WV! i 6 s 
V V —— cos‘ (5 ) 
Fyi" Dys + Ewa (57) 
Fat = Dis! + WVcos$ * W V'cos) - $'W Vein? 
E g E 
"Doy t€" ue x o nec en (58) 
g W V 
à 6 


Consider now the derivatives to be evaluated for the Euler in y 
for the unbounded case; the following appear: (1)Fy,; (2)Fyi, 


(3)Fy From (55) 


F= D+ WV V'cosý + Wsin% 
g 
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Whence: 
Fy - Dy - W V V'sinë t Ncos X (59) 
8 
Ey * Dy: (60) 
Fy! Im Dy! i (61) 


These, then,are the desired derivatives as functions of drag or of 
drag derivatives. The drag derivatives will be determined after a 
look at the derivatives appearing in the bounded equations. From 
equations (43) and (46) for the bounded case we have the following 
derivatives: 
(1) W' 
(2) wi 
(3) fy 
(L) fy 


(5) £' e 
Now W' = Cp fw V! 


gE(f-D-WsinY) ' (62) 
From which 


Ws = CT fW y! fi + Wi + Y os it" - DI - W'sin¥ ay) 
g(f-D-Wsint)} £ W W (£-D-Wsin?) 


(63) 
Since the thrust, f, used in the bounded equations is an actual 
thrust, (i.e. thrust available), engine data for the specific 
engine =» airframe combination will. supply the information necessary 
for evaluating the derivatives on thrust, 
f = f(V,h) A known relationship for a given aircraft 
Then 


fy” of/dV This can be taken directly from engine curves (6ha) 
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f 


Now 


f, = ðf dh (6hb) 


$ 
Sh SY 
and of/dh is available from engine curves while òh% is a simple 


geometric relationship. 


f!» df Vi + df y! 6hc) 
dir. seel 


To evaluate the derivatives on drag it is most convenient.to 
express drag in the following form: 


om D+ By 


= 1 ps VOp, + 1 ps W_ Al 
2 2 sr et 


Or, substituting for Cj its equivalent in terms of lift from equation 


(21a): 
De 1 pë vec, + 2 (WcosY)* (2 SË e ? (65) 
dr Te pS V E | 


In taking derivatives on the above equation for drag, the subsonic 
region only will be considered, No Mach effects will be present 

then and both Cp, and e will be treated as. constants. An extension 
of this analysis to include Mach effects will be found in Appendix B, 


Thus 











Dy = PD» - ] (Weos3)* l} E Y = E 
TeAR pS v3 E 
= 2Do- 2D (66) 
V V 
Dyr E 0 (67) 


Dy+*'=-0 (68) 





3^ 


Dy = -l W cos% a + ei? 


Temes E 
= a? D, tanY (69) 
/ y2 1) 
D: = L (Wcost)2 (1 + — (E) 
«Tel ps ve g 
di (70) 


ye 
Dy! z Ii? cos? y (a + De zur » 23! tanJ $ (veg 4 231 VW!) P 
Ter ps Y JE + 3! p) 7 


a 2D; QW! = 2y'tany = (9' * V (vg" « 2V'4") 
*s') * | (71) 


t€ dë 


DI = Dy V' + Dy $! + Dy Eer + Dy W! : f i 


_ tye 
= 2D yt + t - QD. e 
04 A] i 


= " > Wi + 1 (72) 
way FF 






vi + ' tan) 


The atmosphere parameter, fa , and its derivative, ye ', may 
be evaluated as follows: 

For altitudes below 36,089 feet: 

je Po = ker 1076 m )42% 

/ e 
pit -he256 „(1 - 6,89 x 1078 n )3-?56 (6,89 x 10-Ó dh ) 
A ds 

But dh stand. 

Gr 
Then: 


3 


f 4 
F— = = 29.30 tan Y (73) 
jë 10° 6.89h 
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For altitudes above 36, 059 feet: 
—(h - 36084 


f . pozo) e e 


ui (BROER) x 
pë Lr" g Fever — 


20736 ° 
and 
e. = - tany 
p 20706 (7h) 


/ 
The formulation of the final equations for the unbounded 


case may now be made as follows: 


Fy =- F Tua! -K!u. ( y'tan¥ - y' | 40 
| Ww y 

Substuting for ry, H/V, Fyr', and Fy: from (55), (56), (57), 
and (53); and for Dy weich oceurs in Fy from (66) gives after rrarransing: 


Po - pj (a U 
TR ETE yy ) 
g 


Ed 





“i EO ye 


suite self ut + P + a) val) = 0 





P 


(75) 
Trom (35), the second Buler equations: 


Py + Ptan¥ - Fy! - rg Vea u ) = 0 


V 
Substituting for Fy, Ftant , Fyl! and Fy from (59), (60) ama (61) 
an? for the drag derivatives which are introduced from (69), (70) 


and (71) sives after rearranrins: 


(Dj DRE * 5 - 24 2 ae! 
COSË — A x") M P 
- JV Vu" « VEE) 7 | 
7 U | > 
à (1 + a] (76) 
E 





Solving (75) for V! and (76) for" yields the two equations 


whose solution will define the optimum path, 





"a 
= ge | UE + wor 
WOr V^ cos Y "E =) e 
n Vj tani 
PE - +SS + agaj, ten (77) 
Gje që ct të 
= + y" E B. m Y. 4 " d" a. 
ur en J| 8r 
mo a pi) - ul (78) 
“ve 8 5 —— 


The final equations for the bounded case are formulated in 


like manner: 


From (l3), the first Euler equation 


we on fy - d ve = 0 
Wi | f - Wsin WT 


Substituting for W"/W' from (62) and (63) and for Dy and D' which 


then remain from (66) and (72) yields on rearranging: 
s z£) 
[3 el +2 Dy Ge di 
V V 
l + 
" Me j 


23 tant + y i € Ww 2 - f’ + Wisin? 


+ WA! cosy + | fto V'fy + wW) (f - D - Wein‘) z 0 (79) 





"I E W 


From (46), the second Euler equation: 


EXI HL. y N (f -D- Wsin?) ~ fy - Dy = Weos Y + Dytë 
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‚ (ft! - DI - Wei -WY'cosY) -0 


2p 
i —D-Wsin 


Substituting for W"/W' from (62) and (63)s for Dyi $ Dyi'g and D' 
from (70), (71), and (72) yields on rearranging: 


ft + 3W! « V" — 2(f! . Wisin’ - WY'tcosy) - 2D! 


ae MY Ey - D - Wsın 
pe + x! 

+ Yo. fy (f = D-Wsing) - f, + 2Dj tanë - cos! | 
“etm 1f | 


+ V(VYN + 201 11) - f! -2*'tan] =0 (80) 
GER TË w 


Solving (79) for ¥" and (80) for V" yields two equations 
whose solution will define the desired path: 


ye = i 1 ig - W'siní - W Í'cosY = E ECW Due * i) (£-D-Wsind) 
2D4 f 





T W 


ZUM Fr) "(4 +) ce 
g 


5 


* vet gf £r (81) 


ve 9. Vf! 23w V! 2U' (f! zW'sini - W j'cos)) - 2V!D! 
LIN W - D - Wsin 


J 


s a ce : 3") E (SD Wsind) -- Tet 2D; tand - voed 


- ve (Vin say MO) e vip +v dtan (82) 


is EE 





PART V: METHOD OF SOLUTION OF T48 SQUA TIONS FOR 
AN EXTREMAL PATH 
Considering first the unbounded or general equations for an 
aircraft flying a minimum fuel path, the following equations control 
the aircraft's motion: (21b), (65), (77) and (78). These equations 
will be repeated here for ready references 


A. The general equation of motion along the flight path: 


F=D+ W V! Vcosd + Wsin8 (A) 
E 
E. The drag equation: 
D Ef VG, + _2(Weosd)* (2 + g (B) 
EX) TERE ve E 


C. The first Euler equation: 


Y! os g Do (2 - Yon - Dy 
WO[V^ cos 3 8 


ca 


ETS tan) 
d N nu Or 


D. The second Euler equation: 


i ex)? | * | las 
3 (+ v) | Gem vi ix P) ar 


w n - a | (0) 


£ 2 
Da Aap?’ Cp, 


2 = 2(WeosQ)2 E € oy! 2) 
es‘ 8 





where 


This system is composed of four nonlinear differential equations 








in Y and Ve Three of the equations are first order and one is second 
order. This is equivalent to a third order equation in one unknown 
and thus three initial or boundary conditions may be specified. 

Since no general closed form solution exists to this system, we 
are forced to integrate by means of a step-by-step process. This 
requires the use of initial conditions rather than boundary conditions, 
From equations (A) and (B) it can be seen that the forces on the air- 
plane depend only on V, V', Y , and %' and therefore our initial 
values should be given in terms of three of these at a given altitude, h. 
Since V' can also be determined directly in terms of the other three 
from equation (C), it will Le most convenient to specify V, and 
¥' initially. 

For a given aircraft we know Wọ and Cp,» and for the unlounded 


case we assume unlimited thrust and knowledge of CT at any thrust, 
Now, if we specify V, ¥, and `N! initiaily, V' is computed from (C) 
and D, Dg and D4 from (B). With V' and D known, F is determined 


from (A). Since Cp is a function of F, it was necessary to estimate 


initially to find V'. Now that F has been determined, an iterative 
procedure should be used to correct V', Cr and F. W! is found easily 
from: x 


WY sm =- Cp F (E) 
Vcos 


Now 4" can be computed from equation (D). Thus we have defined all 
quantities for our starting point by initially specifying three 


variables V, y y and Mu. 





By using the eee two terms of a simple Taylor series exo 
pansion we may now integrate stepewise to point number 1 and so on in 
the following manner: 

Knowing Vo', Yo!» Wo' and X o" we compute by Taylor series 

Vy * X4 * V Ax 

Ma = Xo+ Yoh Ax (F) 

t JE! + do Bx 


WL = Wo + Wo! Ax 


" Also since 





ht = tan% (F) 

E = mr Af 
Now with Vj, B1s %]'s W], hy and equations (A) through (E) we can 
determine all quantities at point (1). Using Taylor series again we 
proceed to point (2) and so nn. 

If greater accuracy is required than obtained by two term 
Taylor series, additional terms of the series may be used, the Milne method 
may be used, or the increment in x may be decreased, Decreasing 
the increment in x reduces the errors due to non-linearities of 
the function involved. The Milne method is described in reference 2 

In performing a step-by-step solution to the unLounded case, 
several considerations should be kept in mind to insure that the path 
so determined will be a "flyable" ones 

1.) The equations can possibly ET maximum fuel path instead 
of a minimum. 

2.) Flight path angle, Y „cannot be allowed to drive altitude to 


negative values or the a/c into stall. 





3.) Accelerations and decelerations must be within the physical 
capaLilities of the aircrafte 

4.) Minor oscillations of the computed path in its initial 
Stages are probably due to comoutational variables such as the 
Size of x and to the initial conditions, and can be faired out for 
an actual flight pathe 
Juggling of initial conditions and increment in x will help avoid 
these difficulties at the beginning of the path, 

To compute a minimum fuel path for the bounded case equations 
(C) and (D) must be replaced by (81) and (82). Equations (81) and 
(82) will be restated here only in functional form: 

First Euler Equation 

a, X, Vaa (G) 

Second Euler Equation 

f(V, 375 V', Yt, Ve ADO (H) 

It must be remembered in the bounded case that thrust is 
fixed as a known function of V and Y for a particular aircraft. 
With this in mind, a close imspection of all of the equations for the 
bounded case (i.e. (A), (E), (G), (H)) tells us that the second Euler 
equation (H) is not necessary to the solution of the system, This 
is as we would suspect from the physical conditions of the bounded 
problem, since, by setting thrust at a fixed rating, say maximum 
continuous, it becomes a known function of V and Y. Thus one less 


equation is required for solving the sustem, 





113 


This can be seen more clearly from a look at (G) and (H), 


with the balance equations. We have two second order differential 


equations with second order terms on Y in both (G) and (H) 

equations and second order terms on V in (G) only. y" could be 
Eliminated between (G) and (H) and our system would be equivalent to 
a third order equation in one unknowne Thus three initial conditions 
may be specified. Looking at (G) alone with the balance equations 

we also have an equivalent third order equation in one vaftiable, 

Thus three initial conditions may be specified in either case, 

If we initiall; apecify V, y, andy! as before, for a particular 
aircraft, F, Cp, and D are determined. Then since F and D are known, 
y" y. determined from equation (A)o If V' is determined from 
equation (A), the V" term in equation (H) is redundant. The other 
variables in equation (G) are functions of known or specified 
variables so V, Y, V', y' andy" are related uniquely in (G) and 
therefore (H) can be disregarded. Ey using only the first Euler 
equation all title are specified or determined for the initial 
point of the bounded problem. Fy equations (F) representing the 
first two terms of the Taylor series, ve can again integrate stepwise 
to point 1 and so on as before. 

ı The bounded problem and also the more ganeral unbounded problem 
each have specific applications which fit completely within the scope 
of one or the other of the two situations. It is felt, however, in 
the absense of any computed results that the general solution will 
be^ composi te of the two. Due to the similarity of the systems of 


equations for the two cases, transition from one case to the other is 





uh 


an easy matter at any point on the path. When thrust required ex- 
ceeds that available while using the unLounded equations, transition 
to the bounded system of equations should be made. This transition 
should be accomplished at the Last point on the path at which de 

_ thrust requirement is equal to the thrust available. 

Now, while proceeding rin along the path computed by 
the bounded system of equations, it is necessary that each path point 
‘be checked in the unbounded eqdationse At that point on the path 
where the unbounded thrust requirement comes within the envelopeof 
thrust. available from the power plant, transfer back to the unbounded 
system of equations should be made and computations continued until 
such time as the thrust requirement again forces transition to the 
bounded case. This procedure is continued until the required distance, 


X is covered, 











CONCLUSIONS. 


It is concluded that a system of equations for maximum range 
or minimum fuol flight paths maj Le written which are not restricted to 
level flight or to negligible accolerations. The equations can be 
epnlied to flight paths which commence at take-off or at any other set 
of initial conditions from which an extremal path is desired. 

It is further concluder that L3 specifyinz three initial 
tions, methods of solution bi stenwise intesration exist for 
the sistem of equations. 

t is also concluded thet, shoul1 the extremal path call for 
greater thrust than is availlalle from the aircrait, a maximum range or 
minimum fuel path can be commuted Lased on maximum continuous thrust 
availatle, There are specific arrlications for Loth types of paths 
sadividually as well as for a cornosite patie 

From Appendix A it is concluded that classical parameters for: 
maximum range are restrictive applications of the more general equations 


derived in this thesiso 
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APPENDIX A 


Reduction of Euler Equaticns to Classical Range Parameters 





Froa the Euler equations representing a maximum range path at 
variable altitude it is a simple matter to derive equations or 


parameters for maximum range in level flicht. Starting with equation Ga 
v = ga) PSDo (1 - Or V Y a?cos (i9 V Y'WNwt-3- 
Cr 2cosY g A p V3 P WU V 


E ETC zit) - Wtan?[l + pe! 
g B 








Imposing the conditions for level flight (¥ = ¥' = 0) we have 


4 
remaining: 
bos ad Ee p - =) w h de z-| (4.1) 
Wor 2 g ag de g 


Substituting for W^ its level flight equivalent (ëm pov) 


Bl" 1 psv SL" * 1 Psp 
waa N 2 vel 2 








Then | 
V'WOT »  psCp - 07 V]+ 1psVEny /- 3 - Or 
a es ( 2 zl ri 2) 
And l 
VVO C V x= = + - Cr E uf + Cp V (2:92) 
gë g g 


Assuming norinal values of Op, V, and Cj, say: 
Cp = 11b/lb br = 1/3600 sec 
V = }00 ft/sec 


C = de 


* in slightly different form 


17 


Then Cp V * ,0035 and maj be considered negligible with respect to le 
6 
And V! Cy Cp vé = ,009 V' where V' is in ft/sec and is an extremely 
ft 
8 


Small quantity. This term may be considered to be zero. With these 
k 
assumptions the equation reduces to: 


CD, = 30D, = 3% = 3%, 


Cp = 3 Do (4.3) 


Equation (145) represents a parameter defining the speed schedule for 


v 


t 
maximum range at any one given altitude, This parameter is perhaps not 


as familiar as the classical 
Ci = K p A. 
/ cp 7 l ) 


parameter but it can be shown that the two are equivalent, 


Differentiating cp = K 


dc, velt 
Then 
K = to 
Oy * 
and V ! 
2 
Gi = qef 
Cp ley, me 
from which 
Cp =» 1402 = ho, = U Cp - bo, 
Tek 
and 


Cy = + €», (A.5) 








APPENDIX B 


Consideration of Mach Number Effects 


In the foregoing developement of maximum range and minimum 
fuel equations, only the subesonic Speed region was considered. In 


this section, Mach effects on the equations will be determined and 


' applied, thereby removing any speed limitations. 





Mach number appears explicitly in the drag terms of both the. 
unbounded and bounded cases, and may passibly appear in the thrust, f, 
of the bounded case. Therefore, by determining the Mach corrections 
to the D and f terms and their derivatives, the complete system of 
equations is corrected and is applicable to the entire speed range. 

Considering drag first, the Mach effects thereon may be accounted 


for as follows: 
D = 1 pSVCp 
>f 


= i psvteo + 2(WcosJ)? (a + euy 
2 


ww el pSV g 


Now CH, and e are functions of Mach number. Defining 


Cp, “ vU 
e =: $01) 
Drag may be — asi " I op) 
D = ae ye) + 2(Wcost) = 
2 di see) 


Now 


Do = 1 pS VnM) 
O +f 1 


Di = 2(Wcost)* ( l + yest 2 
1 





Taking derivatives on drag: 


dy mt do. +. Bp Jm 
Y OV "SM "OV 


oD = AD syn + dD ð 
oM ser OF OM 


Since M = V/V, 


am’ = 1 
OV Ve 
Therefore 


2/4 4. M V? _ yy? 
3p. * psvyan » h(Weos)*(1 * 252 fa m 
mr TR sp Y (m) 


Zj 
E i Me | ) IE 
E ST MAR pS v2 TEG) SE Vo 


Or 





Since neither drag nor Mach number are functions of V', as before: 


om 


the derivatives on drag with respect to 1 and 4! are unchanged 


gD = O , also since òM = OM = O 
ng 


from their subsonic form provided CD, and e are inserted as proper 


functions of Mach number. 


To evaluate Dy', fran (70) and the above discussion: 
h(Weos3)*/ 1 + y'vé 
ms (M) 
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From (71) and (B2) 
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To evaluate D'! 
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Then from (pl) and (72) 
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As was previously noted, the bounded thrust, f, may possibly 
be a function of Mach number. If this is the case, derivatives 
may be taken in the normal manner as follows: 


f e f(V, h, M) 
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The partial derivatives with respect to the various variables being 


Obtained from known power plant characteristics. 
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